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Using optimized perturbation theory, we evaluate the effective potential for the massless two 
dimensional Gross-Neveu model at finite temperature and density containing corrections beyond 
the leading large- iV contribution. For large- iV, our results exactly reproduce the well known 1/JV 
leading order results for the critical temperature, chemical potential and tricritical points. For finite 
JV, our critical values are smaller than the ones predicted by the large- JV approximation and seem 
to observe Landau's theorem for phase transitions in one space dimension. New analytical results 
are presented for the tricritical points that include 1/JV corrections. The easiness with which the 
calculations and renormalization are carried out allied to the seemingly convergent optimized results 
displayed, in this particular application, show the robustness of this method and allows us to obtain 
neat analytical expressions for the critical as well as tricritical values beyond the results currently 
^— ; ■ known. 

(N 
•*-> 
O 

O 

t> 

I. INTRODUCTION 

The study of symmetry breaking/restoration in quantum field theories is an important issue of relevance in many 
areas of physics. For example, today, problems regarding phase transitions in Bose-Einstein condensates (BEC) or in 
• Quantum Chromo dynamics (QCD) concentrate a lot of theoretical as well as experimental efforts. 

Of topical importance regarding studies of phase transitions in quantum held theory is the reliability of perturbation 
theory and its eventual breakdown. For instance perturbation theory at high temperatures breaks down due to the 
\^ ' appearance of large infrared divergences, happening for example in massless field theories, like in QCD close to 
critical temperatures (in field theories displaying a second order phase transition or a weakly first order transition 
""^5 ■ 0, 0)j or just because that at high temperatures there are parameter regimes where conventional perturbation 
schemes become unreliable when powers of the coupling constants become surmounted by powers of the temperature. 
In these cases, a nontrivial problem arises since non-perturbative methods must be used. Various non-perturbative 
' techniques have been used to deal with these problems. Among them we will be particularly interested in the 1/JV 
approximation Q which, here, will be considered mainly for comparison. Though a powerful resummation method, 
the 1/JV approximation can quickly become cumbersome after the resummation of the first leading contributions, 
. rH ! like for example in JV = 2 (e.g, BEC and polyacetylene) or JV = 3 (e.g., QCD) finite JV problems. This is due 
to technical difficulties such as the formal resummation of infinite subsets of Feynman graphs and their subsequent 
5h I renormalization. 

An alternative non-perturbative analytical method that we will make use in this work is the optimized perturbation 
theory, or linear S expansion (LDE) 5]. In calculational terms (including renormalization) its appeal regards the 
fact that one remains within the framework of perturbation theory. Then, non-perturbative results are obtained 
by optimizing the perturbatively evaluated quantities. This procedure amounts to eliminate, variationally, mass 
parameters used to deform (interpolate) the original action. Recently, the LDE has been successful in treating 
scalar field theories at finite temperature and/or density. Relativistic scalar Xcj) 4 theories have been treated at finite 
temperature [|| as well as at finite temperature and density 0. At the same time, their nonrelativistic counter part, 
which is relevant for BEC, has been studied in connection with the problem of the dependence with interactions of 
the critical temperature shift AT C , given by the difference between the interacting and ideal gas critical temperatures 
(for a recent review see for instance Q and references therein) . It suffices to say that the latest LDE results are in 
excellent agreement with the Monte Carlo results, considered the most accurate prediction for AT C , and it performs 
much better than the 1/JV expansion used at leading and at next to the leading order (for a review on the different 
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results and methods used to study the AT C in the BEC problem, see [ill). Moreover, the convergence properties for 
this critical theory have been proved by the present and other authors 0, 0, . 

The present work will focus on the LDE applications to problems involving phase transitions in asymptotically free 
models at finite temperatures and densities. Here, we shall consider the 1+1 dimensional Gross-Neveu model (GN), 
which is extensively used as a prototype model in studies related to phase transitions in particle physics as well as 
condensed matter physics. One recent application concerning particle physics is the next-to-leading order in l/N 
evaluation of the effective potential at finite temperature performed in Ref . |l3j . In the condensed matter domain, a 
very interesting recent application of the GN model was to the study of polymers |14| , where a massive version of the 
model was used to show the appearance of a kink-antikink crystal phase that was missed in a previous work M 

In Ref. ^3 the authors work only within the leading order of the l/N expansion, so their results still show a 
phase transition at finite temperature T and density fj,. On the other hand they consider the case of inhomogeneous 
background fields, which should be physically more relevant in the 1+1 dimensional GN model. This is so because, in 
one space dimension, due to the Landau-Mermin-Wagner-Coleman theorem 0, 0, 0] , no phase transition related 
to a discrete symmetry breaking (in this case a discrete chiral symmetry in the massless GN model considered in this 
work) is expected at any finite temperature. In the GN model this can be explained by the role played by kink-like 
inhomogeneous configurations [19( that come to dominate the action functional, instead of just homogeneous, constant 
field configurations. By accounting for kink-like configurations in the large- N approximation, the authors of Ref. 0] 
(see also Ref. [2(J for a review) find evidence for a crystal phase that shows up in the extreme T ~ and large [i 
part of the phase diagram. The other extreme of the phase diagram, large T and /i ~ 0, remains identical to the 
usual large- N results for the critical temperature and tricritical points, which are well known results mmin. in 
the study performed here, we only consider homogeneous backgrounds, but go beyond large- AT, so the phase diagram 
changes as a whole. However, we cannot see any crystal phase at small n, that should be a consequence of kink 
like configurations dominating the action functional at that extreme of the phase diagram. Despite of that, we are 
still improving the calculations for the GN model even though we are not considering inhomogeneous fields. And, as 
explained above, since inhomogeneous backgrounds do not seem to change appreciably the large T and small /i region, 
up to the tricritical point of the phase diagram, we are certainly improving the knowledge in that part of the phase 
diagram. At the same time, wc believe that our results are not faithful in the small T and large fi part of the phase 
diagram, which gets affected at larg e by inhomogeneous backgrounds as shown by the results obtained for the GN 
model in the large- N approximation [141 ], though no results beyond large- A^ are currently available. Eventually, in the 
future it would be opportune to contrast the results found by the authors of Ref. |14j who considered a inhomogeneous 
background field to evaluate the effective action with the ones provided by the LDE in the same context so that it 
could generate the effects of kink configurations beyond the large N limit considered in Ref. . 

Another purpose of the work done here is to show the advantages and reliability of an alternative non-perturbative 
method like the LDE in the understanding of the phase diagram of the massless GN model when considered beyond 
the large- A'' limit. In this case, the massless GN model provides an excellent testing framework for the following 
reasons. As mentioned in the previous paragraph, large-A" results for the critical temperature (T c ) [2l|, critical 
chemical potential (/i c ) |22|. and tricritical points fl5| are well known j2j|. At the same time, Landau's theorem 
|!6lll9| | for phase transitions states that they cannot occur in one space dimension, so that rigorously T c = 0, meaning 
that the large-A" approximation behaves poorly in this case. These two extreme results allow us to gauge the LDE 
performance in connection with the problem since we know that for A" — > oo our results should converge to the 
"exact", although wrong, large- AT result. For finite N, on the other hand, our results should predict smaller values 
for the critical temperatures, in accordance with Landau's theorem. At zero temperature and density, the LDE has 
been applied to this model with some success 0> E3 j since in this simpler case the LDE could even be summed 
to all perturbative orders (at least in the 1/AT approximation), so that the large order behavior of the LDE could 
be inves tiga ted. In fact, its convergence properties have been proved for a particular perturbation series in this 
context [23. At finite temperature, an early application to the GN model [2(| showed the potentiality of this method. 
However, the renormalization program has not been addressed in Ref. [2(|. Here, our aim is to use all the latest LDE 
improvements to evaluate the GN effective potential at finite temperature and density for any value of N. 

This work is organized as follows. In the next section we present the model. In Sec. Ill we review the well 
established l/N results to leading order considering the following four situations: (a) T = fi = 0, (b) T ^ 0, /i = 0, 
(c) T = , fi + 1 and (d) T + 1 , [i 7^ 0. In Sec. IV we present the LDE method and the interpolated GN model, 
evaluating the effective potential for A" — > 00. We show, in accordance with Ref. (2?|], that when correctly applied, 
the LDE exactly reproduces, already at first order, the large- Af results. The situation is unchanged, at any order in 
(5, provided that one stays within the N — > 00 limit. This nice result is valid for any parameter values. In the same 
section, we explicitly evaluate the l/N correction that also appears at the first LDE nontrivial order. In Sec. V the 
LDE order-<5 results are presented and compared to the large-A" results for the four situations described above. Our 
major result is the production of analytical relations for the fermionic mass, critical temperature, critical chemical 
potential as well as tricritical points containing a finite A" correction. We show that all these quantities depend on 
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an optimized mass scale set by the LDE. All the analytical expressions have been cross checked numerically. In Sec. 
VI we contrast the LDE and the 1/N approximation results to leading order and to the next to the leading order. 
Our conclusions are presented in Sec. VII. Two appendices are included to show some technical details and the 
renormalization for the interpolated model. 



II. THE GROSS-NEVEU MODEL 

The Gross-Neveu model is described by the Lagrangian density for a fermion field ipk (ft = 1, ... , N) given by [28 



C = tpk (i ft) i>k + rriFipkipk + y (V'fcV'fc) 2 > ( 2 -l) 

where the summation over flavors is implicit in the above equation, with e.g. ipkipk = SfcLi V'fcV'fc- Since we restrict 
to two-dimensional space-time dimension, ipk represents a two-component Dirac spinor for each value of the flavor 
index ft. When mj? = the theory is invariant under the discrete transformation 



i/j - 75^ , (2.2) 

displaying a discrete chiral symmetry (CS). In addition, Eq. i|2.1|) has a global SU(N) flavor symmetry. 

For the studies of the model Eq. (|2.1I) in the large- TV limit it is convenient to define the four-fermion interaction as 
g 2 N — A. Since g 2 vanishes like 1/N we study the theory in the large-iV limit with fixed A (see e.g. pol^. 

At finite temperature and density, we can study the model Eq. (|2.1|) in terms of the grand partition function given 
by 



Z(fi, /i) = Tr exp [-j3 (H - M Q)] , (2.3) 

where (3 is the inverse of the temperature, /i is the chemical potential, H is the Hamiltonian corresponding to Eq. 
(|2.1() and Q = J dxtykl^ipk is the conserved charge. Transforming Eq. 12.3|l to the form of a path integral in the 
imaginary-time (Euclidean) formalism of finite temperature field theory |30j |. we then have 



Z(0,n) = J nf =1 LW^ fc exp{-S £ fe,V> fc ]} , (2.4) 



where the Euclidean action reads 



SEbPk^k 



dr 1 dx 



ipk {P + M7o - rn F ) ip k ~ Tr^i^kipk" 2 



2AT 



(2.5) 



and the functional integration in Eq. (|2.4I) is performed over the fermion fields satisfying the anti-periodic boundary 
condition in Euclidean time: %1>}.(x,t) — —i))k{x,T + (3). 



III. REVIEW OF BASIC 1/N RESULTS TO LEADING ORDER 



Let us briefly review some of the standard large- N results for the GN model. Considering m F — in Eq. I|2.1II we 
start by looking for a fermionic mass that can be generated dynamically via radiative corrections. This exercise will 
also allow us to set up notation, conventions as well as reviewing useful formulae to be used within the formalism of 
finite temperature and finite density. As usual, it is useful to rewrite Eq. (|2.1|) expressing it in terms of an auxiliary 
(composite) field a, so that 29] 



£ = ~ipk(i ft) ipk ~ cripkipk 



a 2 N 
~2\~ 



(3.1) 
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By using the solution of the equation of motion for a in Eq. (|3.1|l we recover the original model. For renormalization 
purposes, we can also add to Eq. (|3.1[) a counterterm Lagrangian density, £ ctl whose most general form can be 
expressed as 



£ c t = ipk (iA ft) i/) k - BcriPkiPk - C—- + D^ k ip k + Ea + X 



(3.2) 



where A, B, C, D, E and X are renormalization counterterms with the latter representing the zero point energy 
subtraction. 

The appearance of a non- vanishing vacuum expectation value for a, (a) — a c 7^ can be associated to a mass term 
for the fermion field. This is better studied in terms of the effective potential for a, V e g(a c ). As it is well known, 
using the X/N approximation, the large- TV expression for the effective potential is [29| 



V cS , 

N 



2A 



d 2 p 
(2tt)2 



ln(p 2 



(3.3) 



The above equation can be extended to finite temperature T and density fx using the usual associations and replace- 
ments (see Appendix A for details and notation), with the result 

v N a 2 r 

-2&(a e ,T,n) = £-J [ Wj) (a c )+Tln(l + e X p{-[a; p (a c )+/x]/T}) 

+ rin(l + exp{-[ Wp (a c )- M ]/T})] , (3.4) 
where lo 2 (<j c ) = p 2 + a 1 and J p denotes integration over space momentum. The T = 0, fj, = term in Eq. I|3.4|) gives 



V N a 2 r 

-^(a c ,T = 0, M = 0) = ^-J u p {a c ) 



2A 



e< E M z 



2\ e/2 



2A (4tt) 



1 + ln 



(47T) 

M 2 



0(e) 



(3.5) 



where M is the arbitrary mass scale introduced by dimensional regularization. The divergent term in Eq. (|3.5|l can 
either be rendered finite, by imposing a renormalization condition directly on the effective potential by defining a 
renormalized coupling constant \r, via 



da 2 



N 



(3.6) 



or by direct use of the counterterms in Eq. (|3.2|l . in which case we only require a mass counterterm for the a field 
that exactly cancels the divergent term in Eq. I|3.5|l . The two renormalizations are of course equivalent, only differing 
by a different choice of mass scale in l|3.5ll . Choosing the latter form of renormalization, we immediately get the 
renormalized relation 



V N „ a 2 1 



1 + ln 



M 7 



(3.7) 



From the T, [i dependent term of Eq. JH2J, 



V N „ 



[Tin (1 + exp {- [u p (a e ) + /i] /T}) + Tin (1 + cxp {- [u p (a c ) - fj] /T})\ , 



(3.8) 



we can take the limit T — > and perform the momentum integral to obtain the finite result 
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V- 



N 
off 



(<7 c ,T = 0,m^0) = — 0(/x-<7 c ) 

iV Z7T 



(3.9) 



The general expression Eq. (|3.8|l for T ^ and [i ^ 0, however, does not have a close analytical expression, but we 
can express it in terms of a high temperature expansion, in powers of \xjT and cr c /T, in analogous form as it is done 
for the bosonic like temperature dependent momentum integrals |3l|. Defining the function ii(a, 6) by 



h(a,b) 



dx 



In (l + e -^+^- b ) +ln(l + e -^+^+ 



(3.10) 



Eq. (|3.8|l can be written as 



V N xFl 
f (a c ,T, M ) = -T 2 I Jx (a c /T, M /T) 

(27T)2r(2) 



AT 



(3.11) 



We now take cr c /T = and ///T = 6 < 1. Expanding Eq. I|3.10|) in powers of a and b, the result is finite and 

given by [§3 



where £(3) ~ 1.202. Using Eq. (|3.12|) in Eq. I|3.11|l . the high temperature expansion for V c g (cr c , T, ji) becomes 



V 



N 



ttT\ a 2 clE , a c 2 7C(3) 



where we have dropped terms that do not depend on a c . Let us now review four important situations 



(3.13) 



A. The T = n — case: the fermionic mass at large-iV 
At T — fi — the large- N effective potential is given by Eq. (|3.7|) , 



V N rr 2 1 

off (~ rp n ,, s ! c ^-2 

(0>, J = 0, Ll = U) = er. 



1+ln 



(3.14) 



We can already check here that, contrary to the classical (tree- level) potential where the minimum occurs at a c = 0, 
one now has symmetry breaking that is quantum generated since the minimum of the effective potential occurs at a 
non- vanishing value given by 



dV£(a c ,T =0,n = 0) 



= 



(3.15) 



where a c sets the large- TV result for the fermionic mass, at T = fj, = 0, as 



?tif(0) = (T c = A/exp ( — — 



(3.16) 
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B. The T/0 and fj, — case: the critical temperature T c 

Physically, the case fi — means that the number of fermions and anti-fermions are the same. For this case, from 
Eq. I|3.7[) and using the result (|3.13[1 for the finite temperature term for the effective potential, we have that 



7C(3) 4 



2tt ' 327r 3 T 21 c ' 

which shows chiral symmetry restoration at a critical temperature T c given by 



(3.17) 



T c = m F (0) ~ 0.567 m F (0), 



(3.18) 



in accordance with Ref. |2J|. This transition is second order, as can be easily checked by writing Eq. (|3.17(l in the 
form 



|K,T # o,,=. ) = i,„(l).; + ||.; 1 

which for T close to T c gives a temperature dependent vacuum expectation for ct c {T) = mp(T) in the form 



(3.19) 



m F {T) 



8tt 2 lT r - T 



7C(3)- 



(3.20) 



that shows a continuous transition at the critical point T = T c . For T > T c the fermion mass mp(T) vanishes, 
restoring the chiral symmetry. 

C. The T — and jj, ^ case: the critical density fi c 

For the case of zero temperature but unequal number of fermions and anti-fermions, // ^ 0, one considers Eqs. I|3.7|l 
and H3.9[) . which give the expression for the effective potential at finite density but zero temperature. Then, 



, M 2 
1 + In — 



o- 2 ln 



(3.21) 



For /i < a c the minimum of the potential reproduces the same result as Eq. (|3.16|) . For fi > a c we have from Eq. 
(|3~2H that 



a c (jjb) = yJm F (0) [2/j, - 771f(0)] , 
which is valid for m.p(0)/2 < [i < mp(0). In the case a c = m_p(0) and from Eq. (|3.21() we obtain 



(3.22) 



V N 1 
^(a c , T = 0, M < <t c ) = ~^m F (0) 



(3.23) 



For \x < ttif(0)/2, V^(ct c ,T = 0,/i ^ 0) will have a minimum for a c = mi?(0) and a maximum for a c = 0. For 
mp(0)/2 < n < mp(0), V^(a c , T — 0, /i ^ 0) has minima at tr c = and a c — mp(0) and a maximum at cr c (/i). There 
is a value of /i where both minima satisfy 



V e %(a c = m F (0), T = 0, /x c ) = ^(^ = 0, T = 0, Mc ) , 



Nt 



(3.24) 



where 
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m F (0) 

Mc = —jf - ' ( 3 - 25 ) 

in agreement with Refs. [22,|2i|]. This result indicates a first order phase transition: for values of /i > fj, c , a c — TOi?(0) 
is a local minimum (false vacuum), while <f c = is a global minimum (true vacuum). For \x < /i c the two vacua get 
interchanged. 



D. The T/0 and ^ / case 



This situation can be analyzed numerically only. On the (T, fi) plane one can identify a tricritical point P tc where the 
first and second order transition lines meet. In units of m_p(0), the numerical result is P tc — (T tc , fi tc ) — (0.318, 0.608) 
[l5| . The large- N phase diagram will be presented in subsection V.D together with the LDE result. 



IV. THE EFFECTIVE POTENTIAL FOR THE INTERPOLATED THEORY 



Let us now turn our attention to the implementation of the LDE procedure within the GN model. According to 
the usual LDE interpolation prescription |£| (for a long, but far from complete list of references on the method, see 
[23) the deformed original four fermion theory displaying CS reads 

Ce(ip, tp) = tp k (i ft) V-fc + v(l - 5)i) k i> k + 5^{^ki>k? • (4.1) 

So, that at 8 = we have a theory of free fermions. Now, the introduction of an auxiliary scalar field a can be 
achieved by adding the quadratic term, 



SN ( A 1 J 



— U+-M k ) , (4.2) 



to £s(tp,tp)- We are then lead to the interpolated model 



Cs = (i ft) ipk - Saipk^k - vi 1 - S)ipkipk - T^e 2 + £ct,s , (4-3) 

where the counterterm Lagrangian density, C c t t s, has the same poly nomial form as in the original theory, Eq. (|3.2|l . 
while the coefficients are allowed to be S and r\ dependent [fil l24| . Details of the renormalization process for the 
interpolated model are given in the appendix. Note that the same interpolation of the form l|4.3|l was also used in 
Ref. 34], but it is different from the ones used in Refs. 0> HE E3- HH • I n those references the interpolation was 
not carried out in the original four fermion theory but on its bosonized version. However, we argue that the present 
choice is more adequate because at 5 = one has only free fermions. Otherwise, the quadratic bosonic term <r 2 /(2A) 
survives at 6 = and the theory looks to be composed by free fermions and bosons but this is misleading since, by 
the equations of motion, a = (X/N^kipk- From the Lagrangian density in the interpolated form, Eq. (|4.3I) . we can 
immediately read the corresponding new Feynman rules in Minkowski space. Each Yukawa vertex carries a factor 
—iS while the (free) a propagator is now —iX/(NS). The LDE dressed fermion propagator is 



p - ry* + le 



where 77* =77 — (77 — o~ c )5. 
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A. The LDE Effective Potential in the Large-A Limit 



In the following we start our analysis of the LDE effective potential in the large- TV limit, i.e., the LDE is applied 
directly to the effective potential after the large-A limit was taken. In practice, as already emphasized, within the 
LDE one perturbs in powers of S only, without having to notice the actual powers of N as we do here. However, in 
this subsection, we do the 6 expansion within the l/N leading order only to be able to compare both methods. In 
other words, we first want to establish the reliability of the LDE-PMS by considering N — > oo in order to further 
compare our results with the "exact" ones furnished by large- N approximation. The reader should not be confused 
by this apparent mixed type of expansion. Using the new Feynman rules from the interpolated model 14. 31 1. in the 
large- N limit, one can write the effective potential, analogous 1 to Eq. (|3.3() . to arbitrary orders in 6 as 



Y^(a c , V *,N^oc) = 6^ + ij -^tThx^-v*) , (4.5) 

which can be seen as obtained directly from the series containing all perturbative terms, Eq. (|3.3|) . simply by 
performing the replacement a c — > 77*. Developing Eq. I|4.5II to first order in 6 one obtains 



V if xi , „cri? f d d p , , f d d p ri~o r , 



N v " ' 2A J {2n) d ^ " J {2n) d i>-r\ 

Fixing 6 = 1, we now optimize V e R for rj as in most of the previous references on the LDE method by using the 
principle of minimal sensitivity (PMS). In the PMS procedure one requires that a physical quantity & k \ that is 
calculated perturbatively to some fc-th order in 6, be evaluated at the point where it is less sensitive to this parameter. 
This criterion then translates into the variational relation 361 



d$(fc) 
dr/ 



#7,5=1 



. 



(4.7) 



The optimum value fj that satisfies Eq. I|4.7|l must be a function of the original parameters, including the couplings, 
thus generating non-perturbative results. In the following we will apply Eq. (|4.7|l directly to the effective potential. 
Applying the PMS to Eq. I|4.6I) immediately gives the result fj = a c> which then recovers exactly the large- N result, 
Eq. I|3.3I) . This same trend holds at any temperature and/or density, value of the coupling A, as well as number 
of space-time dimensions since there is no need to perform the integrals explicitly. Also, as shown in Ref. [2]j, the 
inclusion of higher order terms does not spoil this nice result, provided we stay within the large- N limit, since they 
are all of the form 6 k fk(il — cr c ) k , where k > 2. Note that the PMS admits another solution given by 



d 
dr] 



d d p 



tr- 



r) + ie 



= 



(4.8) 



However, this solution, which depends only on scales introduced by the regularization process, and thus is not 
proportional to the basic scale of the model after dimensional transmutation: A/e _7r / A , can be taken as unphysical 
|27j since it brings no information about the theory being studied. 



o 





+ 




+ . . . 



In this subsection we shall work in an arbitrary number of space-time dimensions to emphasize the generality of our optimization 
procedure. 
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FIG. 1: Diagrammatic relation between Eqs. 14,51 and 14. 6H . The thick continuous line on the LHS represents the LDE dressed 
fermionic propagator which is rj* dependent. The first diagram on the RHS is an order-<5° vacuum graph. The dashed line 
represents the auxiliary field, a c , while the black dot stands for the 5r/ vertex. The last two graphs on the RHS are of order-5. 



One can spot a subtle point associated to the LDE evaluation of the effective potential by looking at the diagrams 
considered to order-<5 in Eq. 1)4.6) 1 . These graphs are displayed in Fig. ^ which shows the diagrammatic relation 
between Eqs. (|4.5|) and 1)4.6)1 , Note that Eq. 1)4.6) 1 contains two a c independent (vacuum) terms that would be 
neglected in most evaluations since they are irrelevant as far as CS breaking/restoration are concerned. However, as 
pointed out in Ref. [27^ . they are r\ dependent an so must be considered until the theory is optimized. One can also 
easily see that, in fact, those contributions are responsible for the quick LDE convergence towards the exact large- N 
result already at the first non-trivial order. The explicit form of the LDE large- N effective potential is quickly obtained 
by applying the finite T and fi rules described in the previous section to Eq. (|4.5() . In particular, from Eqs. (|3.7|l and 
1)3.11)1 for d = 2, one then obtains the renormalized result 



00)=^-- 



1 , M 
- + In — 

2 V 7 ?* 



(4.9) 



Expanding Eq. (|4.9|) to order S at T = \i = one retrieves, as before, the PMS result fj = cr c that reproduces the 
large- N result. The same optimum value for 77 appears when we express Eq. I|4.9|l in the other limits studied in the 
previous section. At the same time, as discussed above, the PMS gives another solution, Eq. 1)4.8)1 . which is fj = M/e. 
Rather curiously, this result coincides with a c = mp(0) evaluated at A = it. 



B. The Effective Potential to Order 5: going beyond the large- TV limit 



In the following we show how the optimization procedure implemented by the LDE improves over the large-iV 
results. Here we revert to the usual LDE procedure by expanding the effective potential in powers of 5 only. This 
quantity can be expressed in terms of the full fermionic self-energy, whose terms contributing up to order S 2 are 
displayed in Fig. [3 The corresponding contributions to the effective potential up to order S 2 are shown in Fig. 
[3J This nicely illustrates how the LDE incorporates, at the same perturbative order, graphs that have different N 
dependence. 




FIG. 2: Contributions to the fermion self-energy up to order S 2 . The continuous thick line represents the LDE, rj t dependent, 
dressed fermion propagator while the thin lines represent the r\ dependent propagators. The dashed line represents the scalar 
auxiliary field, a. Tadpole diagrams are not shown since they do not contribute to V e s 
. The second graphs represents a correction to the a propagator while the third has corrections to the Yukawa vertex. The 

first and fourth are exchange (rainbow) type of graphs. 




FIG. 3: Feynman diagrams contributing to the effective potential at order 5 2 . Note that the first and second have a 1/7V 
dependence while the last two have a 1/N 2 dependence. The thick lines in the first two graphs represent the LDE, 77*, 
propagators which, as discussed in the text must be further expanded to order-5 2 . 

Let us now consider the second graph of Fig. [3]which, for comparison purposes, has been neglected in the previous 
subsection. It reads 
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V 



(a) 

:ff,gl 

N 



(v* 



d 2 p 

W) 2 



tr 



1 — 77* + ie 



(4.10) 



where the trace is over Dirac's matrices only 2 . The term E represents the first contribution of Fig. [3to the fermion 
self-energy, 



£ Q (?7*) 



d 2 q 



(2n) 2 



it 



(4.11) 



Originally, Root 37] has summed up all 1/N corrections to next to the leading order and by further perturbatively 
expanding his result one easily retrieves Eq. i|4.10|) . 

Then, using the rules described previously, at finite temperature and density, one finds, after taking the trace and 
renormalizing with the fermionic mass counterterm (see the appendix), the result for Eq. I|4.10|l 



v: 



(«) 



cff.5 1 , rf, v x A 



N 



Air 2 N 



rj* 



+ T 2 J 2 ( V */T,fx/T) 



(4.12) 



where 



I 2 (a,b) = -2 



dh(a,b) 
da 2 
dx 



\/x 2 +< 



1 



,\/x 2 +a 2 +b _|_ ^ e ^/x 2 +a 2 — b _|_ ^ 



with Ii(a, b) defined by Eq. H3.10|) where, now, a = f}*/T, and similarly 



(4.13) 



J^ipbjb) = sinh(fe) / dx 

JO <<)S 



h(va? 



cosh(fo) 



(4.14) 



The (finite) contributions given by the term proportional to l2(a,b) and J2{a,b) originate from the summation over 
Matsubara frequencies of the T and /i dependent contributions, more precisely from Eqs. i|A2|l and (IA3|) respectively 
(see Appendix A for details). Note that the divergence is only contained, at this order, in the T = 0, fx = part, 
which is renormalized by standard counterterms in the MS scheme (see Appendix B for details). 

Finally, by summing the contribution Eq. (|4.12|) to the effective potential expression to leading order in N, Eq. 
(14. 9|) , and expanding the resulting expression one obtains the complete LDE expression to order-<5 
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In 



+ S 
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4tt 2 N 



M 
1 

M 



2T 2 I 1 ( V /T,fx/T) 
I 2 ( V /T,»/T) 



ln( — )-I 2 {n/T,n/T) 



+ T 2 J 2 ( V /T,v/T) 



(4.15) 



Notice once more, from Eq. I|4.15|l . that our first order takes into account the first next to leading order correction to 
the large-A^ result. 

When considering the case T — 0, [i ^ in Eq. (|4.15|) one can take the limit T — > 0, which then gives for the 
functions I\, I2 and J2 the results 



2 The factor —1 corresponding to a closed fermionic has already been taken into account l37l . 
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Urn T a I 1 ( v /T, t i/T) 
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o(m - v) 



2 / m + vV - 2 

77 In - //yM _ »r 



Hm7 2 (77/T, M /T) = %-77)ln 



lim TJ 2 (r]/T,n/T) = sgn(//)0(/j - 77) ^M 2 - V 2 



(4.16) 

(4.17) 
(4.18) 



Note that all the above results vanish for fi < rj. When considering the case T ^ and /1 = the high temperature 
limit, t)/T <C 1 and fi/T <C 1 will prove to be useful. For Ji, this approximation follows from Eq. Ij3.12jl . while for I2 
it can be obtained using Eq. H4.13JI which yields 



I 2 (a, b) = In (1) - 7E + (V + y) + 0(« 4 , & 4 ) 

In the case T ^ and /i ^ the integrals 2i, I 2 and J2 will be handled numerically. 



(4.19) 



V. OPTIMIZATION AND NUMERICAL RESULTS BEYOND LARGE- TV 

Before proceeding to the specific d = 2 case, considered in this work, let us apply the PMS to the most general 
order-<5 effective potential which is given by 



— -7T- (c c , v) — S— r + i f , n ?. tr In (i> — 7?) + Si [ , n f , tr- — — 

N 2\ J (27r) d KF " J {2ir) d p-rj + ie 
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(2Tr) d p 1 — rj + ie 



(5-1) 



This exercise will help the reader to visualize the way the LDE-PMS resumms the perturbative series. After taking 
the traces in Eq. I|5.1|l one obtains 
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Then, setting S — 1 and applying the PMS one gets 
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(5.2) 
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Po \ d ( . f d d p 



(2ir) d (p 2 - if + ie) J dr\ 
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l + ri— 
dr\ 

Po 



d d p 



1 



(2ir) d {p 2 -r/ 2 + ie) 



{2n) d {p 2 - t? 2 + ie) 



(5.3) 



As one can see in Appendix A (Eq. (|A3|l ) the last term of the equation above only survives when /i ^ 0. In the case 
fi = 0, Eq. (|5.3|l factorizes in a nice way which allows us to understand the way the LDE-PMS procedure resums the 
series producing non-perturbative results. With this aim one can easily check that (at S = 1) 



£ (/X = 0,T),T) = -j=T) 

Then, when /x = the PMS equation factorizes to 



d d p 



(27r) d {p 2 -r, 2 + ie) 



(5.4) 



[77 - a c - £ a (/j = 0, 77, T)] ( 1 + 77— 



(27r) d (p 2 - t; 2 + ie) 



(5.5) 
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which leads to the self consistent relation 



77 = <7 C + E a (n = 0,7?, T) 



(5.6) 



which is valid for any temperature and number of space-time dimensions. In this way the LDE fermionic loops get 
dressed by a c as well as rainbow (exchange) type of self energy terms like the first graph of figure|3| Typical optimized 
Feynman graphs are shown in Fig. The graphs in part (a) of that figure represent the order-5 contributions prior to 
optimization while the infinite set shown in part (b) represents the non perturbative optimized result. Note how the 
graphs are dressed by rainbow type of self energy contributions. This was expected since at order-i5 the perturbative 
LDE effective potential receives information about this type of topology only. If one proceeds to order-5 2 information 
about corrections to the scalar propagator as well as Yukawa vertex (see figure |2J) will enter the perturbative effective 
potential. Then the PMS will dress up these perturbative contributions and so on. In other words, the simple 
evaluation of a first topologically distinct graph will bring non perturbative information concerning that type of 
contribution. Figure 0] clearly shows that the LDE-PMS re-sums all powers of 1/N corresponding to the (rainbow) 
class of graph. 







(a) 






(b) 

FIG. 4: (a) Feynman diagrams contributing to the effective potential at order S prior to optimization, (b) The infinite set of 
graphs contributing to the optimized effective potential. 



Note that the mathematical possibility 



d d p 



(2ir) d (p 2 - fj 2 + ie) 



= 



(5.7) 



corresponds to the unphysical solution discussed in Sec. IVA (Eq. I|4.8[) h 

Having illustrated the way the LDE-PMS resummation works, let us concentrate in the d = 2 case by collecting 
our results for the complete order-<5 effective potential, Eq. (|4.15|) . once the optimization equation is applied to it. 
Using the PMS procedure we then obtain, from Eq. (|4.15|l at 77 = fj, the general result 



XT 2 d 
+ 2^V J2(? ' /T ' /i/T) ^ J2(? ' /T ' /i/T) 



77=77 



0, (5. 



where we have defined the function 



yu,.T.,n = In ( M) - I.ui/T.t'/T) 



(5.9) 
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Eq. H5.8|) expresses our general PMS result, Eq. (|5.3() . for the specific d — 2 case. This can be easily seem by recalling 
that in this number of space dimensions (and 6 = 1) the exchange (Fock) type of self energy is given by 



E (?7, T, fi) = -V^j^yiv, T, m) 



(5.10) 



In the following we will present the analytical (whenever possible) and numerical results for Eq. I|5.8|l . For 
convenience, the results will be presented in units of M for different values of A and N. We start by analyzing the 
simplest case, of zero temperature and density. 



A. The T = fx = case 

Taking Eq. (£>3 at T = /i = [that is, ^(77, T = 0, fi = 0) = ln(M/r])] one gets 
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1 



77 In ( — 







(5.11) 



As discussed previously, the first factor leads to the coupling independent result, 77 = M/e, which we shall neglect. 
At the same time the second factor in (|5.11H leads to a self-consistent gap equation for 77, given by 



%!(o-c) = Cc 
The solution for 77,51 obtained from Eq. 1)5. 12[) is 
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2ttN \M 



(5.12) 
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2ttA ct c 
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exp 



2ttAT\ 



(5.13) 



where W(a;) is the Lambert (implicit) function (23, which satisfies W(x) exp[W(x)] = x. 

To analyze CS breaking we then replace 77 by Eq. (|5.13|l in Eq. I|4.15|l . which is taken at T = and fi = 0. As 
usual, CS breaking appears when the effective potential displays minima at some particular value a c 7^ which is 
obtained by minimizing the resulting effective potential with respect to a c , 



9V eSS i (a c , 77 = rjgi 



8=1, 



— + - 77 In — = . 

A 7T M 



(5.14) 



Since mp = cr C) after some algebraic manipulation of Eq. Ij5.14|l using the definition of the W{x) function, we find 



m FS x (0) = a c (T = 0, fi = 0) = M^(A, AT) 1 



2A^ 



(5.15) 



where we have defined the quantity J-(X, N) as 



J-(X, N) = exp 



A[l - 1/(2A)] 



(5.16) 



Eq. H5.15|l is our result for the fermionic mass at first order in 6 which goes beyond the large- N result, Ea. (|3.1ti|l . 
Note also that in the N — > 00 limit, since N — * 00) = exp(— tt/A), Eq. (|5.15() correctly reproduces, within the 
LDE non-perturbative resummation, the large- TV result. In Fig. Elwe compare the order-<5 LDE-PMS results for a c 
with the one provided by the large- A'' approximation. 

At the same time, 77 evaluated at a c (T = 0,/i = 0), to first order in 8, also follows analytically from Eqs. (|5.13|) 
and H5.15[l . 



rf S1 (a c ) = (1- — ) a c 
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(5.17) 
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FIG. 5: The dimensionless minimum <r c (in units of M) as a function of A for T = fi = 0. The dashed line represents the 
TV — > oo result while the continuous lines were produced by the LDE-PMS at order-5. The numbers beside the curves identify 
the value of N for each case. 



The analytical results Eqs. (|5.15|l and (|5.17J) are exact expressions (of course, at this first order of the LDE) following 
only from the minimization condition Eq. I|5.14|l after some algebra. In particular, we emphasize the simple scaling 
relation obtained between fj and a c in Eq. (|5.17|) . that only depends on N, leading to the result (|5.15|l . and that 
will prove to be a very good approximation even for the more general case, T ^ 0, as we shall see next. Note that 
this relation is the explicit form, for this case T = ji = 0, of the more general relation Eq. I|5.t)[) above, namely: 
Z a ( f i = 0,r 1 ,T = 0) = -a c /(2N). 

Actually, there is an alternative simpler way of deriving the results (|5.15fl - <|5.17l) without using the solution for fj in 
Eq. inHSJ, that will moreover prove useful later on when we shall consider the more complicate situations with n =/= 0, 
or typically also if we would consider higher LDE order contributions. That is, instead of using Eq. I|5.14|l as giving 
a c (fj), we can use it to substitute the rather complicated logarithmic dependence \nfj/M in Eq. I|5T2|I . canceling out 
at the same time the A dependence, and obtaining thus a very simple linear equation for a c /fj which only depends on 
N, 



H 1 + *T) ■ 

which immediately leads to the results H5.15fl - I|5.17|) . The solution in Eq. (|5.13(l is useful however, as it gives fj for 
any a c values, away from the minimum in Eq. i|5.14H . 



B. The T 0, ^ = Case 



For the finite temperature but zero chemical potential case, still using Eq. I]5.8p . the optimized rj is now determined 
by the solution of the expression 



y( V , T, f i = 0)+rj—y(7 ] ,T,fi = 0) 



i] - a c + rj 



2ttN 
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(5.19) 



The solution coming from the first term in Eq. I|5.19|l . corresponds to the unphysical optimized result, whereas the 
solution obtained from the second term gives the equivalent of the self-consistent gap equation, (|5.12|) . 
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-l 

(5.20) 

We could next proceed by numerically solving Eq. I|5.2Q(I for fj and substituting it in place of r] in Eq. I|4.15|l (evaluated 
at fi = 0). Then, by minimizing the effective potential we would obtain the general behavior of a c as a function of the 
temperature. The critical temperature for CS restoration would then be determined by the solution of a c (T = T c ) = 
as usual. However, here an explicit analytical result for T c can also be obtained if we apply the high temperature 
approximation to Eq. (|4.15l) . with r//T <C 1, and then optimize the resulting expression. The validity of using the 
high-T approximation before the optimization procedure could be questioned in principle since 77, at the level of Eq. 
(|4.15l) . is arbitrary. However, one may easily perform a cross check by performing a numerical PMS application, 
as described above, without using the high-T expansion. The results we found for T c from both approaches agree 
very well with each other, showing that the high-T expansion is valid in a large range, though actually rj and T c are 
numerically of the same order of magnitude. A simple reason for this is that the true high temperature expansion 
parameter is not r] 2 /T 2 , but rather 7 /A[((3)/(2tt) 2 }j] 2 /T 2 ~ 0.05?7 2 /T 2 , as it is clear from Eq. (|53T|) below, so that 
even if T c ~ rj the expansion parameter is still small enough. If we then expand Eq. (|4.15|) at high temperatures, up 
to order r] 2 /T 2 , we obtain 
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Then, setting (5 = 1 and applying the PMS one gets 
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(5.22) 



where, once more, the first term in Eq. I|5.22|l corresponds to the unphysical (A independent) solution, while the 
second term gives 



ri(a c ,T) = a c < 1 + 
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(5.23) 



whose solution can be expressed in the form 
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(5.24) 



The result given by Eq. (|5.24|) is then plugged back into Eq. (|5.21|) for 5 = 1 and the resulting expression is again 
expanded for <j c /T <C 1 in the high-T approximation. The order of the transition can easily be checked numerically 
simply by plotting the effective potential for different values of T as shown in Fig. 

The result shows that the first order LDE result for finite N predicts a continuous phase transition as in the large- N 
case. By extremizing the effective potential, at high-T, we obtain a maximum at a c = and two minima at: 
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(5.25) 



Figure shows a c /M given by Eq. (|5.25|l as a function of T/M, again showing a continuous (second order) phase 
transition for CS breaking/restoration. 
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FIG. 6: The effective potential, V e fi/N as a function of a c for /i = and T = 0.170757 M. The parameter values are A = 7r and 
TV = 3. The continuous curve, which shows CS restoration, was plotted using the LDE optimized results. The dashed curve, 
which signals CS breaking, corresponds to the large- TV predictions. Both V e s/N and a c are in units of M. 



0~c 




0.15 0.17 0.19 0.21 

FIG. 7: The non-trivial minimum a c as a function of T for the parameter values A = n and TV = 3. The first order LDE curve 
(continuous line) displays a continuous phase transition occurring at the critical temperature T c — 0.170 M while the large- TV 
result is T c = 0.208 M. Both a c and T are in units of M. 



The numerical results illustrated by Figs. HJ1 and [7\ show that the transition is of the second kind and an analytical 
equation for the critical temperature can be obtained by requiring that the minima vanish at T c . From Eq. (|5.25|) 
one sees that a c (T = T c ) = can lead to two possible solutions for T c . The one coming from 
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(5.26) 



can easily be seen as not been able to reproduce the known large- TV result, when TV — * oo, T c = M exp(7^ — n/X)/ir. 
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However, the other possible solution coming from 



f Me lE \ 

-2A^+(27V_i)Aln(— —\ 
jives for the critical temperature, evaluated at first order in S, the result 
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The 



with T\(N) as given before, by Eq. (|5.16|) . Thus, Eq. (|5.28|l exactly reproduces the large- N result for N — » 
results given by this equation are plotted in Fig. ISJas a function of A for different values of N. 

As already mentioned, the critical temperature predicted by Eq. I|5.28[l coincides, to a very good approximation, 
with the one observed in the numerical results starting from Eqs. H4.15fl an d l|5.19|) without the use of the high-T 
approximation. We note in Eq. I|5.28[) the very same scaling relation than the one for tof(O) at zero temperature, in 
Eq. I|5.17[) . thus involving a relation that only depends on N . This simple scaling derives in fact from the equation 
for the optimized solution fj(T), where to obtain the factorized form Eq. I|5.22|l . we have neglected C(l/T 4 ) small 
terms. Plugging back these neglected terms, and manipulating Eq. H5.22J) . one easily obtains the relation: 



0-c = v(T) 



(5.29) 



where we defined for convenience the high-temperature expansion parameter x = 7/4[C(3)/(27r) 2 ]^ 2 /T 2 ~ 0.05 r? 2 /T 2 , 
so that x 2 (\/Tr)fjH (fj, N) defines the remnant part of order 1/T 4 , whose explicit form we do not need to specify. 
Thus, neglecting these (D(x 2 ) terms in Eq. I|5.29[) gives the very same relation between fj and a c than at T = in 
Eq. (|5.17(l . and furthermore, Eq. (|5.29|) is directly related, after some algebra, to the result H5.28J) for T c above. 
Accordingly, though the simple scaling relation is generally not expected to hold at arbitrary T ^ 0, the deviation 
from this relation is essentially negligible due to the smallness of those 0(x 2 ) corrections. We shall argue later on, in 
section V.E, that this universal scaling property is expected to remain a good approximation as well at higher orders 
of the LDE approximation. 

The (non-perturbative) LDE result shows that T c is always smaller (for the realistic finite N case) than the value 
predicted by the large- N approximation. In the light of Landau's theorem for phase transitions in one space dimen- 
sions, which predicts T c — 0, our LDE results, including the first 1/N corrections, seem to converge to the right 
direction. 



C. The T = 0, fi Case 

As discussed in Ref. 0] this extremum of the phase diagram is the very sensitive to the role played by kink 
like configurations. However, in the present work only homogeneous background fields are considered and we are 
not in position to compare our results for this part of the phase diagram with the ones provided in that reference. 
Nevertheless, we are in position to contribute by computing finite iV corrections so that one may, eventually, use the 
LDE-PMS in conjunction with inhomogeneous background fields to further improve the phase diagram found by the 
authors of Ref. hj. The case of zero temperature but finite chemical potential (density) also follows from Eq. (|5.8|1 . 
Using Eqs. (|4.16() and l|4.17[) in Eq. (|5.8() . we find two situations. In the first, for 77 > /.i, the optimized fj is found 
from the solution of 

For rj < /i, using the relations (|4.1t)|) , 14.17fl and l|4.18|l in Eq. (|4.15|) we obtain that the optimized ?y is the solution of 
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FIG. 8: The dimensionless critical temperature, T c (in units of M), as a function of A for T / and n = 0. The dashed line 
represents the N — > oo result while the continuous lines were produced by the LDE-PMS at order-5. The numbers next to the 
curves identify the value of N for each case. 



The solution for Eq. (|5.30|> is exactly the one obtained previously, given by Eq. (|5.13(l . Concerning Eq. I|5.31[) one can 
again obtain an analytical solution, by following a reasoning similar to the one done for T = fi = leading directly to 
Eq. H5.18f) . but with slightly more involved algebra in the present case. Thus, consider the non-trivial minimum a c , 
obtained from dV c g/da c ~ 0, now for /i ^ 0. From Eq. H4.15[) it follows that a c is given by 
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(5.32) 



which replaces Eq. I|5.14|l for /i ^ 0. Now again we can use this to eliminate simply the complicated logarithmic 
dependence ln[(/z+ vV 2 — V 2 )/M] in Eq. (|5.31|l . thus obtaining after straightforward algebra a second order equation 
for a c as a function of fj and /i, whose explicit solution reads: 
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(5.33) 



where 



G( „,,, W)s i( 1 __L) (5 .34, 

contains the fi dependence. The relation (|5.33|) is the appropriate generalization, for /x ^ 0, of the simple scaling 
relation obtained at T = fi = in Eq. I|5.17|l . We have eliminated the other possible solution for a c (namely with 
yCTT — > — in Eq. (|5.33|) '). by noting that for A — > 0, Eq. I|5.33|l correctly reproduces the simpler scaling relation 
in l|5.17|l . while the other solution would give a c — > 0. 3 It will prove also useful to expand l|5.33|l in powers of A: 



3 The leading order large N relation u c = r\ is also consistently reproduced for TV — > oo in Eq. 15.331 . 
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which should be valid thus for moderate values of A/V. This immediately shows that the corrections due to fi ^ to 
the simple scaling obtained previously for T = fi = in Eq. I|5.17|l are actually suppressed by O (A 2 / (tt 2 N)) , moreover 
the \i dependence enters only at the next order, (\/ir) 3 . These properties are somewhat analogous to the case of 
the T (/i = 0) corrections to the scaling in Eq. (|5.29[) which are also suppressed by the small high-temperature 
expansion parameter. We shall come back in the next sub-sections to the important consequences of these relations 
for the general case T,[i ^ 0. 

Next, we can extract the critical chemical potential, /i c , in an analogous way as in the large-TV problem shown in 
Sub-sec. IIII CI and, as in the large- N problem, we can check that a first order transition also occurs. To this aim we 
first calculate the effective potential at the value a c given by Eq. (|5.15() . where the relevant expression is given by 
^eff(cr c ) for /i < rj. After some algebra many terms cancel out so that we obtain simply 



V 

1 - l/(2N) 



V£\* c = * c ,r, = fj) = -^, (5.36) 

i.e., the same expression as the leading order N — > oo one, but with the appropriate fermion mass at first order in 
S, fjgi given by Eq. i|5.17|l . This has now to be compared with the value of V c s f° r ^7^0 but with 77 = 0, which is 
simply obtained from Eq. (|4.15|l as 



Vtf' (<* = 1 = 0, /,) = -£ (1-5^) , (5-37) 
so that we can deduce an analytic expression for the critical density /i c , 



valid at this first order in S. The appearance of a pole in Eq. (|5.38|) for A = 2irN is an artifact of our first order in 5 
approximation and thus probably not physically relevant. The point is that going to higher orders in the expansion in 
S it will also bring different corrections of same order in powers of 1/N, if expanded. Thus, we could for example limit 
ourselves in this analysis to the first 1/N order in this expansion, since complete 1/N 2 corrections are not included 
at the order in S we are considering. Seen this way, we could also expand Eq. (|5.38l) in powers of 1/N, 



M 
71 



exp(— 7r/A) 



tt A 
2AA + AttN 



0{1/N 2 ) 



(5.39) 



which then exhibits no pole. 

We also proceeded numerically to obtain the solution from (|5.31(l and subsequent values of /^ c , in order to have a 
useful crosscheck of the more complicated most general case with both T and /i finite, where analytical expressions 
are not available. An example for fixed values of N and A is presented in Fig. |5J which shows the effective potential 
as a function of a c . The analogous condition to Eq. H3.24|l leads to the result fi^s 1 — 0.232 M, in the first order of 
the LDE, that should be compared with the large-A result, fi c ,N-nx> — 0.260 M. 



D. The /i/0 Case 



We now arrive at the point of analyzing the complete LDE phase diagram, at order-<5. In the lack of possible fully 
analytical solutions in this general case, we use numerical routines to determine both curves of second order and 
first order phase transitions and their point of intersection which then defines the tricritical point. In Fig. 1101 we 
show the large- TV result compared to the LDE first order result. It can be seen that in the LDE non-perturbative 
approach the region for CSB is diminished in an appreciable way. In units of M, the LDE result for the tricritical 
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V ef f /N 
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-0.007 
-0.008 
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FIG. 9: Large- A (dashed line) and first order LDE (continuous line) results for the effective potential, V e s/N. The parameter 
values are N = 3, A = ir and T = 0. The effective potential has been evaluated at the (first order) LDE critical value 
/i c = 0.232 M for which the large- A approximation still predicts CSB. Both V c g/N and /i are in units of M. 




FIG. 10: Large- A (outmost line) and LDE (innermost line) phase diagrams. The parameter values are A — 3 and A = n. The 
dots represent the tricritical points. The large- A result for the tricritical point is Pt c ,N->oo ~ (T tc ,fJ,tc) = (0.117,0.224) while 
the LDE approximation gives P tc s i ~ (Tt c , Ht^s 1 — (0.091, 0.192). In both cases, the lines above the tricritical point represent 
second order phase transitions while the ones situated below represent first order transitions. All quantities are given in units 
of M. 



point, with N = 3 and A = n, is P tc s 1 = {Tt c , Mtc)* 1 = (0.091,0.192) while the large-iV approximation gives 
PtcJV-oo = ( T tc Mte) = (0.117, 0.224). 

It is worth remarking how the LDE tricritical point falls, approximately, over a line joining the large- N result and 
the origin. This remarkable result is shown in Fig. ^]that shows Pt c for A = n/2, tt, 2ir and N = 3, 10, and N — > oo 
(LN). A deviation of about 9% is observed only for large values of the ratio X/N (e.g., 2n/3). For ratios close to the 
unity the deviation is very small (about 3%). 
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FIG. 11: The tricritical points Pt = (Tie, fitc), in units of M, for different values of A and N. The stars represent A = 7r/2, the 
dots A = 7r and the diamonds A = 2n. The corresponding values of N are directly specified by the numbers on the graph. The 
angle 7 represents the quantity T tc //x tc as discussed in the text. 



E. Generality of Scaling for the LDE Results 

For moderate values of X/N, we thus note the almost invariance, up to very small corrections in the LDE finite N 
case, of the ratio T tc j \itc that defines the angle 7, 



tan 7 = — ~ 0.523 . (5.40) 
(He 

All our results show that this ratio remains to a good approximation largely independent of both N and A. As 
already discussed, this strongly suggests to postulate a very simple (but approximate) result for the predictions of 
the tricritical points and for all other dimensionful quantities obtained for the GN model within the LDE. The result 
(|5.40l) indicates that T tc and fitc, as a function of N and A, to a good approximation in the relevant range, scale 
according to 



T tc (X,N)~cg(X,N)M , (5.41) 
Mtc (A,iV)~ 5 (A,iV)M, (5.42) 

where c = tan 7 and g(\,N) is a function of the parameters A and N. At the first LDE order, we have obtained 
explicitly g(X, N) — J 7 (X, N), where T(N) is given by Eq. I|5.16|l . up to very small corrections in the T ^ (fi = 0) 
case. For the case /x 7^ (T = 0) we saw that the approximation is valid only for A < n approximately, which is easily 
understood by examining the approximate expanded form Eq. 15.35fl , of the relation Eq. I|5.33[l . As anticipated, one 
can see that the corrections for fi 7^ to the simple scaling relation l|5.17|l . strictly valid for T = fi = only, are 
moderate as long as A is not too large, which remain essentially true also for the more general case T ^ and fi ^ 0. 

Interestingly, the actual result for the tricritical points, if compared with the large- N result shown in Subsec. IIII Dl 
is of the form 



Ptc = {Ttc, fhc) =s (0.318, 0.608) fj s t (a c ) , (5.43) 

where fj$i — J 7 (X, N)M . The numerical deviations from the above result are very small (less than 5%) for X/N < 1.3. 

The result (|5.43|l together with our previous results obtained within the first order in the LDE, Eqs. Q5.15[l . (|5.28|) 
and (|5.38|) . respectively for the fermion mass (<7 C ), critical temperature (T c ) and critical chemical potential (/i c ), show 
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the same approximate scaling as given in terms of the optimized LDE quantity 7751 . Also, as we have already seen 
previously, in the large- A limit, 77,51 (<r c ) — > mp(0) and all our LDE results correctly reproduce those of Sec. 11111 

The reason why most of our results (except perhaps for the somewhat extreme case T = 0, /i 7^ 0) exhibit ap- 
proximately universal scaling properties as a function of fj, to a very good approximation, are well understood at 
the first LDE order, as explained in the previous subsections. On more general grounds, in writing the interpolated 
Lagrangian density Eq. (|4.1J) . a explicitly dimensionful (mass) parameter, 77, is introduced in the originally massless 
model where there is no other dimensionful parameter. Thus, from simple dimensional analysis, clearly all physical 
quantities should scale with 77, or actually fj derived from the optimization procedure. 4 More precisely, at first LDE 
order, and T = /1 = 0, we obtained the exact simple relation Eq. (|5.17(l between fj and a c , which only depends on N. 
More generally, in terms of the basic scale of the model, say A = Me~ n ^ x in the MS renormalization scheme, truly 
non-perturbative results are expected to give for the ratio of the relevant quantities mp/A, T c /A and /i c /A, some 
specific dimensionless coefficients depending only on N. A priori, there is no reason why these coefficients should 
be exactly the same for the three quantities. In contrast, on purely perturbative grounds, for T / and fi =/= 
we expect at finite LDE order, to obtain after optimization more complicated scaling relations, i.e., not only with 
different coefficients for the three relevant quantities, but with such coefficients being some non trivial (dimensionless) 
functions of T, fi, <j c , A, and N. Indeed, an explicit example that this is the case already at first LDE order is 
illustrated by the relations Eq. 1)5.29)1 for T ^ 0, and Eq. (|5.33(l for /j ^ 0. Both relations strictly depend on A, and 
T (or, respectively: /_*). However, we have seen that, quite remarkably, this extra dependence upon the coupling is 
quite suppressed, in such a way that the fermion mass mp, critical temperature T c , and even /i c have to a very good 
approximation identical scaling factors which only depend on N. Moreover, we have shown that in the T ^ (/i = 0) 
case, this result is not a numerical accident but can be well understood analytically by noting that the corrections to 
this universal scaling are intrinsically quite negligible. Even in the most general case when both T ^ and fi =/= 0, 
the results for the tricritical point happens in a regime where the simple scaling relation appears to be still a good 
approximation. The only exception is the somewhat extreme case when T = 0, fi =/=: in this case, if A is sufficiently 
large, the appropriate generalization given analytically by Eq. 1)5.33)1 . leads to substantial deviations from the simple 
scaling relation for T = \i = in Eq. 1)5.17)1 . 

We believe in fact that this trend generalizes as well to arbitrary higher orders of the LDE expansion, as we shall 
argue heuristically next. To examine what is happening at higher LDE orders, let us first consider the T = [i = 
case. It is simpler to generalize the reasoning used to obtain Eq. ()5.18|) at first S order, i.e. using the solution a c 
of Eq. 1 )5.14) 1 to replace the logarithmic dependence hi(fj/M) into the (optimization) equation defining fj, that thus 
generalizes Eq. 1)5.12)1 at higher orders. In this way the latter will be a simpler, purely algebraic equation, for the 
(only) parameter cr c /fj, in terms of the remaining parameters, the coupling A and N. More precisely, at higher LDE 
orders, from general arguments the effective potential will have additional perturbative terms of the generic form 



fc>2 V 7 

with leading, next-to-leading, etc different powers of hi(r)*/M), and cll, etc, the corresponding coefficients to be 
calculated from the relevant Feynman graphs. Next, performing the LDE to some given order and taking the limit 
5 — ► 1 etc, the resulting dependence upon a c and 77 is much more involved than at first order, so that one obtains 
from the CS breaking condition a non-linear relation for a c , generalizing Eq. ()5.17)) . and as well for the optimized fj 
solution, Eq. 1)5.12)1 . Nevertheless, it is clear that these two relations can be used e.g. to eliminate the logarithmic 
dependence hifj/M, thus obtaining a polynomial equation for a c /fj, which only depends on N and A. Actually, if we 
were able to re-sum this LDE series to all orders, we would certainly expect that the only dependence on the coupling 
A in all physical quantities would be entirely included in terms of the basic dimensionful scale of the model e.g. in the 
MS scheme: Aj^g = Me~ n ' x . However, at finite LDE orders, the optimization generally defines a rather complicated 
dependence on A for fj. But we have checked explicitly at the next (5 2 ) order that this induces relatively small 
deviations from the first order scaling relation in Eq. 1)5.17)1 (upon assuming obviously that the unknown coefficients 
Cj appearing in Eq. 1)5.44)1 take generic values of 0(1) ). Next, we saw that the T ^ result gives, at this LDE order, 



c NLL {N)\n k ^ 



(5.44) 



4 Note that though at first the interpolation procedure done in Eq. 14.11 seems to explicitly break the discrete chiral symmetry of the 
model, rj is initially an undetermined parameter. The introduction of the auxiliary field a just makes the optimized r) a function 
of the background scalar field cr c . After optimization fj(cr c ) ~ cr c and so, chiral symmetry breaking and restoration have the same 
interpretation as in the original (non-interpolated) GN model. It is a non-vanishing value for the minimum of the a field effective 
potential, (cr c ) = cr c ^ that still signals the breakdown of chiral symmetry, and r/(<r c ) ~ <r c provides the scale for all physical quantities 
derived. 
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a very small correction to the scaling relation existing for T = 0. For fi ^ 0, a deviation to this simple scaling is 
predicted by Eq. (|5.33|) . but it remains, however, numerically moderate up to relatively large values of the coupling 
A. Now, since as we just examined the T = scaling properties are expected to generalize at higher LDE orders, and 
the LDE generally converges quite rapidly, we can expect that the higher order modifications for T, fi ^ to these 
first LDE order scaling properties should remain small corrections. 



VI. COMPARING THE LDE AND THE l/N RESULTS 

Let us compare, in this section, the LDE leading order results with the ones given by the 1/iV approximation at 
leading order (LO) as well as next to the leading order (NLO). As already pointed out, the effective potential, at 
T = and /i = 0, for the 2d Gross-Neveu model to the next to leading order (NLO) was first evaluated by Root 37]. 
The NLO correction to the fermionic mass, at T = /x = 0, was explicitly evaluated by Forgacs, Niedermayer and Weisz 
|39|. Using a combination of the thermodynamic Bethe ansatz and the l/N effective potential at T = fj, = 0, Chodos 
and Minakata ^(j were able to obtain the NLO correction for \i c . Recently, the authors in Ref. ^| computed the 
complete NLO in the l/N expansion for the effective potential at T ^ and /i = performing a detailed numerical 
analysis of their results. They also exhibit in details a number of non-trivial properties, in particular for the expected 
behavior at high temperature, but because of the appearance of a Landau pole near the T ~ T c regime, they do not 
conclude on a well-defined value of T c from the full l/N calculation. It is thus difficult to compare our numerical 
LDE estimates of T c with their numerical results. 

Next, we emphasize that, as far as we know, there are no l/N NLO results for the case T ^ and fi^O, thus no 
results for the tricritical points beyond the large- N approximation. For comparison purposes, let us use our notations 
and conventions to present the only two available analytical results. 

For the case T = fi = Ref. gives 5 




7i 2n 3n 



5 In our comparison with the exact l/N result of the fermion mass as given in Ref. I.Sfll . we should be cautious to remove from their 
expression a term: Mexp(— ir/\)/N, due our scheme being such that our reference scale is M exp(— n/X), rather than the full A^jg- 
expression of I'ill . 



24 



FIG. 12: The fermionic mass mir(0) = a c , in units of M, plotted as a function of A for N = 3, T = and /j, = 0. The dashed 
line represents the 1/N result at leading order, the dot-dashed line represents the 1/N result at next to leading order and the 
continuous line is the first order LDE result. 

For the case T = 0, // ^ Ref. gjj gives 




2 2 

FIG. 13: The critical chemical potential fi c in units of M, plotted as a function of A for N = 3 and T = 0. The dashed 
line represents the 1/7V result at leading order, the dot-dashed line represents the 1/N result at next to leading order and the 
continuous line is the first order LDE result. 

Figure IbH shows the LDE critical chemical potential at T = 0, fi c = M/y2 !F(X, N), as a function of A for N = 3. 
The same figure shows the 1/N results at LO, given by Eq. (|3.25l ), and at NLO, given by Eq. I|6.2|l . As we have 
emphasized in the text, the LDE results for this case are valid only up to A ~ ir. Apart from that, this is the case 
where kink like configurations start to play a major role so any quantitative results must be interpreted with the due 
care Q. 

VII. CONCLUSIONS 

The analytical non-perturbative technique known as the LDE has been applied to the two dimensional Gross-Neveu 
model effective potential at finite temperature and chemical potential. Following the prescription suggested in Ref. 
[27j we have shown that, within the large- N limit, the LDE exactly agrees with 1/N LO "exact" result for any values 
of T or /i. Having established this reliability we have considered the first finite N correction that already appears at 
the first non-trivial order. The variational optimization procedure has produced interesting results that turned out to 
be possible to be cast into analytical form. A careful analysis of our numerical and analytical optimized results has led 
us to write down five simple relations, that take into account finite N corrections, concerning the T = = [i fermionic 
mass (tr c ), [i c (at T = 0), T c (at fi = 0), and the tricritical points, (T tjC , /i t c ). As we have discussed, four of these 
quantities essentially scale, to a very good approximation, with the LDE optimized mass scale at T = = fi given by 
fjgi/M = !F(X,N), with !F(X,N) defined by Eq. I|5.16|l . with fj being the only dimensionful quantity present in the 
interpolated Lagrangian density. The only exception is the case of /i c for T = where, for large enough coupling A, we 
obtained substantial deviation from the simple scaling, with an explicit expression of its dependence on A. However, as 
already mentioned, we believe our results are not faithful in the small T and large \i part of the phase diagram, which 
is the most substantially affected by inhomogeneous backgrounds as shown at leading 1/N order in Ref. [lj|. The 
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basic reasons for the approximately very good validity of the universal scaling relations, at the first LDE order, were 
identified and understood, as explained in some detail in Sec. V. Moreover, we argued that this quasi-universality of 
scaling is expected to hold also at arbitrary higher orders of the LDE. At T ^ and /j ^ our main results concern 
the evaluation of the phase diagram, containing finite iV corrections that, as far as we kn ow, has not been carried out 
before. Comparing our pcrturbative type of evaluation with the ones performed in Refs. 13, 39, 40], for example, one 
may notice some of the LDE advantages. Namely, it automatically introduces an infrared cutoff that makes possible 
completely perturbative evaluations. At the same time, at each (perturbative) order one has just a few Feynman 
graphs to evaluate as compared to the traditional non-perturbative methods, such as the 1/N approximation. In 
particular, this advantage of the LDE procedure means that the renormalization program can be easily implemented. 
Now to proceed beyond the first LDE order, as far as the T ^ 0, /i = case is concerned, it could be interesting in 
principle to exploit those full 1 /N results of Ref. [l3| , by re-expanding these in ordinary perturbation theory in the 
coupling A, and then proceeding with the LDE procedure as outlined in section III. However, after definite efforts to do 
this, it proves to be of very little use for our purposes, simply because beyond the first LDE order one can obtain these 
results only in purely numerical form. While, to perform the LDE, we necessarily need at least to have the analytical 
dependence upon the coupling A and fermion mass to be able to perform the interpolation as described in Eq. 1)4.1 [I . 
Moreover, even at first order, it is very difficult to compare our results with theirs, due to the different momenta 
routing used, resulting in their much more involved expressions of the T-dependent contributions in particular (note 
however that this alternative routing is the only possible one beyond the first order). More generally, at the first LDE 
order here investigated, we cannot expect to get very close to the exact 1/N results, but at the same time we expect 
the LDE to converge faster and it will incorporates terms beyond the 1/A?" at higher orders, as discussed before. 

Speculating on the expected behavior at higher LDE orders, the occurrence of a Landau pole in the complete 1/N 
results found by the authors of Ref. |l3l | therefore invalidating a priori a unambiguous determination of a T c value, 
deserves some general comments. Accordingly, it is useful to recall how the Landau pole emerges in the construction 
of Ref. |13| : in very rough terms, at finite T a pole can occur in the dressed propagator of the a field for some value 
of the T-dependent effective fermion mass, mp{T), because the latter decreases as T increases from zero, reaching 
eventually itif(T c ) = at critical temperature 6 . 

Now, it is interesting to note that the LDE will avoid this Landau pole problem, in a rather trivial way: since 
by construction the LDE stops at finite orders, where there cannot be a pole, while a Landau pole is relevant only 
when considering a resummed perturbative series. Moreover, let us assume that one would manage to re-sum the 
LDE perturbative series to all orders, in some approximation (which can be done explicitly e.g. for the GN mass 
gap |24l l25j at order-l/A, using the LDE together with renormalization group properties): one would thus obtain a 
resummed expression exhibiting possibly a pole at a critical mass value, but even in such a case, the optimization 
prescription used in the LDE construction will escape this pole, i.e the optimization "freezes out" the mass at a value 
which generally cannot coincide with a possible Landau pole value. 

Although convergence properties cannot be accessed in the first non-trivial order, the fact that our T c is, in 
accordance with Landau's theorem, smaller than the LO large N result gives further support to the method. However, 
a deeper discussion about convergence is beyond the scope of the present work due to the technical difficulties 
in evaluating, at T and /j ^ 0, some of the three loop graphs shown in Fig. [3] We recall that the LDE-PMS 
convergence properties in critical theories have been analyzed, by the present authors, in connection with homogeneous 
Bose-Einstein condensates |9J. We intend to extend the present work by considering the order-5 2 three loop graphs 
that will bring the first 1/A 2 corrections. This will help us to gauge convergence properties, the effects of the 1/N 2 
terms, as well as the eventual generalization of the universal scaling relations between the critical and tricritical 
quantities with the LDE optimized mass scale. 
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Not e in contrast that this Landau pole is harmless at T = 0, because in the corresponding exact 1/N calculation of the mass gap (see 
e.g. I3lll4ll ), this pole actually cannot be reached for any consistent value of the mass gap. 
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APPENDIX A: SUMMING MATSUBARA FREQUENCIES AND RELATED FORMULAS 

The integrals encountered in Feynman's graphs are performed, as usual, at finite temperature and density with the 
substitution rules po — > i(u> n — ifi) where fi is the chemical potential and u) n = (2n + l)nT, n = 0,±1,±2, . . ., are 
the Matsubara frequencies for fermions. We sum over the Matsubara frequencies with standard contour integration 
techniques, and regularize the remaining (Euclidean) momentum integrals with dimensional regularization and carry 
the renormalization in the MS scheme. E.g., momentum integrals of functions f{po, |p|) are replaced by (see Ref. 

M) 

For the divergent, zero temperature contributions, dimensional regularization is carried out in dimensions d = 1 — 2e 
and in the MS scheme, in which case the momentum integrals are written as 



dp 



47T 



d d P 

(2^p ' 



where M is an arbitrary mass scale and je — 0.5772 is the Euler-Mascheroni constant. 

The Matsubara sums which are relevant for the different integrals considered in section IV can be derived as (see 
e.g. 123): 
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APPENDIX B: RENORMALIZATION OF THE LDE EFFECTIVE POTENTIAL 

We give here some details on the renormalization procedure for the effective potential of the interpolated model 
(|4.3|) . First, let us consider the non-renormalized result for the effective potential in the large- N approximation, which 
from Eq. I|4.5|l is given by 
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Going beyond large- TV one must add the order-<5 term, Eq. (|4.1U|) . 
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where S a is the term 3c in Fig. O 
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FIG. 14: Feynman diagrams contributing to the effective potential up to order S. Diagrams la, lb, lc and Id are the LDE ultra 
violet divergent contributions. Diagram Id has a symmetry factor 1/2. Diagram 2a represents the counterterm corresponding 
to the zero point energy subtraction. Diagram 2b is constructed with the counterterm (3d) which renormalizes the exchange 
self energy (3c). The linear counterterm 2c can be obtained by renormalizing (with 3b) the direct self energy term, 3a. 



It is worth noting that the relative simplicity of our final result, Eq. (|4.12|) . for this two- loop graph is largely due 
to the appropriate choice of momenta routing, allowing to factorize the result into (squared) one- loop contributions. 
We checked explicitly that it is strictly equivalent to another momenta routing choice in the literature |l3l l37| , more 
appropriate when considering higher order contributions but which at this two-loop level would result into much more 
involved intermediate expressions. 

As it has been shown |(| 0| , the standard MS renormalization procedure and the LDE commute with each other, 
so that one may perform the LDE before renormalization, introducing thus extra (5-dependent counterterms, or 
alternatively directly on renormalized quantities. We shall rather follow here the first approach which may be more 
illustrative for our purpose. An explicit evaluation of Eqs. i|B2|l and (|B3|) with the rules given in appendix A gives 
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Eqs. fBTJ and $B4& 

give the total contributions from diagrams la,lb,lc and Id in figure ITU 
To evaluate the contribution given by graph 2b in Fig. 1141 we need to define the mass counterterm D s exc ipip used 
to renormalize the exchange self energy given by 3c. Note that the divergences are only associated to the terms 
proportional to the mass. This sets the Feynman rule to iD 8 exc where 
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Then, one can evaluate the one loop contribution corresponding to 2b in Fig. [21 
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obtaining 
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so that when adding the two contributions i|B7|l and l|B4|) the 1/e divergence as well as some finite terms cancel out, and 
there only remains the 1/e 2 divergence to be renormalized by zero point (vacuum energy) subtraction counterterms. 

The Feynman renormalization coefficient corresponding to the linear counterterm E s a is obtained from the divergent 
part of the fermionic loop contained in the graph 3a of Fig. 1141 namely 
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Then, adding all contributions one has 
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The remaining divergent contributions come from purely vacuum (field independent) graphs and can be absorbed by 
the zero point energy subtraction counterterm 
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finally giving the finite effective potential 
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which is the result shown in Eq. (|4.15|) . More details about renormalization within the LDE can be found in the first 
work of Ref. (see also [HEI for the GN model). 
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